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The End of the Rainbow: What Can We Say About the 
Extragalactic Sub- Megahertz Radio Sky? 
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ABSTRACT 

The Galactic disc is opaque to radio waves from extragalactic sources with frequencies 
v less than ~ 3 MHz. However, radio waves with kHz, Hz, and even lower frequencies 
may propagate through the intergalactic medium (IGM) . I argue that the presence of 
these waves can be inferred by using the Universe as our detector. I discuss possible 
sub-MHz sources and set new non-trivial upper limits on the energy density of sub- 
MHz radio waves in galaxy clusters and the average cosmic background. Limits based 
on five effects are considered: (1) changes in the expansion of the Universe from the 
radiation energy density (2) heating of the IGM by free-free absorption; (3) radiation 
pressure squeezing of IGM clouds by external radio waves; (4) synchrotron heating of 
electrons in clusters; and (5) Inverse Compton upscattering of sub-MHz radio photons. 
Any sub-MHz background must have an energy density much smaller than the CMB at 
frequencies below 1 MHz. The free-free absorption bounds from the Lyman-a forest are 
potentially the strongest, but are highly dependent on the properties of sub-MHz radio 
scattering in the IGM. I estimate an upper limit of 6 x 10 4 Lq Mpc -3 for the emissivity 
within Lyman-a forest clouds in the frequency range 5 — 200 Hz. The sub-MHz energy 
density in the Coma cluster is constrained to be less than ~ 10 -15 ergs cm -3 . At 
present, none of the limits is strong enough to rule out a maximal Tf, = 10 12 K sub- 
MHz synchrotron background, but other sources may be constrained with a better 
knowledge of sub-MHz radio propagation in the IGM. 



Key words: radio continuum: general - radiation mechanisms: non-thermal 
galactic medium - miscellaneous 
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1 INTRODUCTION 

Whenever a new wavelength window has been opened on 
the electromagnetic spectrum, it has led to new discover- 
ies (|Harwitlll98ll :lL awrencd 120071 ). In the past century, al- 
most the entire electromagnetic spectrum has been explored, 
with few gaps between 10 MHz and 100 TeV. The cosmic 
electromagnetic backgrounds at most of these frequencies 
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1 Although the extragalactic Extreme Ultraviolet (EUV; 100 — 
912 A) radiation is easily absorbed by neutral hydrogen, EUV 
has been detected from AGNs and gal axy clusters at wavelengths 
of < 160 A (see the EUV review by iBowver, Drake. fe Vennesi 
2000). Closer to the Ly man limit, not even a G alactic background 
has been detected yet llEdelstein et al.ll200lT ). and the Galactic 
neutral hydrogen would much more effectively block incoming 
EUV radiation. However, the extragalactic ionizing background 
at a variety of redshifts is indirectly meas ured by the ionization 
state of Lyman-a fores t clouds (e.g., [Baitrik. Duncan, fc Ostrikerl 
ll98Sl: IShuII et alj|l999h . 



have now been measured (as rev i ewed by iRessell fc Turner! 



2006). 



Searches 



Il99d : iFukugita fc Peebles! 12004 iTrimbL 

are underway f or very-high energy photons, inc l uding 
those with PeV (IChantell et alj 1 19971: iBorione et all 11998 : 
ISchatz fc KASCADE Collaboration! 120031 ') and even EeV 
energies (e.g., lAbraham et all l2008bh . These searches 
face several challenges. At PeV energies, extragalactic 
searches are hampered by 77 opacity, in which PeV 
photons interact with the CMB to produce e + e~ pairs, 
and the Universe is highly opaq ue at these energies 
|Moskalenko. Porter, fc Strondl2006|). EeV photons are pre- 
dicte d to exist (e.g.. Iwdowczvk, Tkaczyk, fc Wolfendalel 
Il972l : iGelmini. Kalashev. fc Semikozl 120081 ). and the Uni- 
verse i s possibly transparent to sev eral Mpc at these en- 
ergies (|Protheroe fc Bierm ann 1996), but very low number 
statistics are a problem. Finally, above 10 24 eV, the Galaxy 
becomes completely opaque, as single p hotons pair pr oduce 
e + e~ off the Galactic mag netic field (Stecker 2003). The 
only other frontier in terms of photon energy is at the other 



© 0000 RAS 



2 Lacki 



end of the electromagnetic spectrum, the lowest frequency 
radio waves. 

The Universe is filled with tenuous plasma, the inter- 
galactic medium (IGM), which prevents radio propagation 
below the plasma frequency up — ^/n e e 2 /(7rm e ). Waves be- 
low this frequency evanesce within one wavelength, reflecting 
off the medium. The mean baryonic density of the Universe 
is (rib) ~ 2.5 x 10 _7 (1 + z) 3 cm -3 , although regions in the 
IGM may have higher or lower densities. For IGM with elec- 
tron density 5 e {n,b), the plasma frequency is 

vp = 4.55 e 1/2 (l + z) 3/2 Hz. (1) 

The density distribution in the Universe can be approxi- 
mated as a lognormal distribution, with most of the vol- 

j II 1 

ume being relatively empty ( | Coles fc J ones 1991). Roughly 

90% of the Universe's volume at z = is predicted to have 
5 e > 0.002 (|Bi fc Davidsenl Il997l ). which corresponds to a 
plasma frequency of up = 0.2 Hz. This is essentially the low 
end of the cosmic electromagnetic spectrum, below which no 
electromagnetic wave can ever travel0 

Although there might be electromagnetic radiation with 
u < Hz in the Universe, direct observations at these fre- 
quencies are impossible. Obser vations down to 1 MHz will 
soon be routine with LOFAR (|Rottgerinel l2003) . However, 
the Earth's ionosphere has a typical plasma frequency of 
10 MHz, depending on the time of day and other conditions. 
Similarly, the Moon may have an ionosphere with a plasma 
frequency of a few hundred kHz, which sets a lower limit to 
the f requency of lunar-based observatories (|Jester fc Falckei 
2009). Satellites have measured the Galactic radio emission 
down to ~ 100 kHz (|Brownlll973l ; iNovaco fc Brownlll978l ). 
Space-based observatories near t he Earth, such as the pro- 
posed ALFA (jjones et al.ll2000al lbh. can potentially observe 
down to 30 kHz, which is the plasma frequency of the So- 
lar wind at Earth's orbit. Further out from the Sun, the 
plasma frequency continues to drop, and the Vo yager probes 
took advantage of this to detect kHz emission (|Kurth et al.l 
1 1984 ). However, the interstellar medium itself has a plasma 
frequency of 2 kHz, which serves as a hard limit for direct 
observations of Galactic and extragalactic sources. 

The situation is even worse for extragalactic and distant 
galactic sources. The warm ionized medium (WIM) forms a 
disc with scale height h ~ 1 kpc and a typical density of 
nwiM ~ 0.01 cm. The WIM is opaque to low frequency ra- 
dio emission because of free-free absorption. Although this 
absorption has uses for t omography of the Galactic ISM 
^Peterson fc Webberll2002l ) , it prevents all direct extragalac- 
tic observations at u<3 MHz. There may be a few low 
density 'c himneys' that allow in some lower frequency radio 
emission (I Jester fc Falckell2009l ). but for practical purposes, 
the extragalactic sky at the lowest radio frequencies will be 
shrouded from our direct view for the foreseeable future. 

But there are ways around this limit to locate or place 



2 There is a loophole, since electromagnetic waves below the 
plasma frequency will evanesce with a scale of about one wave- 
length: if the electromagnetic wave has a cosmologically long 
wavelength, it can stretch over the entire Universe. Hawking ra- 
diation from the accelerating expansion of the Universe actually 
does have cosmological wavelengths (e.g., [Gib bons & Hawking] 
1977), but the energy density of Hawking radiation is negligible. 



bounds on extragalactic sub-MHz radio sources. Free-free 
absorption and other opacity sources become more effective 
at low frequency; this is why they prevent sub-MHz radi- 
ation from reaching Earth. Yet, this also means that the 
lowest frequency radio waves are tightly coupled with the 
matter in the Universe. Observations of intergalactic matter 
therefore constrain these radio waves. Synchrotron absorp- 
tion and Inverse Compton scattering also place bounds on 
sub-MHz radio in regions with cosmic rays and magnetic 
fields, like galaxy clusters. I will argue that sub-MHz radio 
emission does not need to reach us, because the Universe is 
our detector. 

I will first discuss our expectations for the sub-MHz sky 
(§ El, including postulated sub-MHz radio sources (§ 12.11) . 
the IGM phases that can interact with sub-MHz radio waves 
(§ 12.20 . and sub-MHz radio propagation through the IGM 
(§[2]3}. I then set new limits on the extragalactic background 
at sub-MHz frequencies. The first limit I consider is the weak 
bound from the expansion history of the Universe (§|3}. The 
second bound is from the heating of the IGM by free-free 
absorption of sub-MHz radio waves (§ HJ. A third bound 
on incident radiation at the lowest frequencies comes from 
the radiation force exerted on IGM clouds (§ [5}. Two more 
bounds can be set for clusters with cosmic rays: bounds on 
synchrotron heating by low frequency radio waves (§ [6} and 
bounds on Inverse Compton upscattered radio waves (§ [7J). 
Finally, if extragalactic photons with energy > 10 20 eV are 
ever detected, they will set extremely strong limits on the 
kHz to MHz radio background (§|BJ. 

My goal throughout this paper is to derive upper 
bounds on the sub-MHz emission with order of magnitude 
accuracy, where possible with current knowledge. In some 
cases even this is not possible - the radio scattering prop- 
erties of the IGM over large distances are not well known, 
and this is a huge source of uncertainty in arguments that 
rely on radiative transfer. Future theoretical work may re- 
duce these uncertainties. Throughout this work, I consider 
the time-averaged sub-MHz background at z £» 0, though 
these methods may be applied to other redshifts. 



2 THE EXTRAGALACTIC SUB-MHZ SKY 
2.1 Sources 

2.1.1 Synchrotron emission 

The MHz to GHz radio emission of galaxies and some galaxy 
clusters is dominated by synchrotro n emission fro m cos- 
mic ray (CR) electrons and positrons l|Condonll 19921 reviews 
the synchrotron ra dio emission of non-AGN galaxies, while 
iFerrari et ail 120081 reviews the synchrotron radio emission 
of galaxy clusters). This emission should continue down 
to lower frequencies, where mostly lower energy electrons 
and positrons are radiating. Non-relativistic electrons and 
positrons emit at the cyclotron frequency 

-=2^= 2 - 8Hz G&) < 2 > 

At GHz wavelengths, the radio spectrum falls off steeply 
as u~ a , where a ~ 0.7 for galaxies and 1 < a < 1.5 for 
clusters, implying greater spectral densities at low frequency 
and greater power at low frequency for clusters. 
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However, the radio spectra of galaxies are unlikely to 
continue down to the sub-MHz range for a number of rea- 
sons. First, CR spectra tend to be flattened at low ener- 
gies, because of escape, ionization, or bremsstrahlung losses, 
or simply because the low energy CR electrons are not old 
enough to be cooled by synchrotron emission. Second, free- 
free absorption will prevent sub-MHz radiation from escap- 
ing galaxies, just as it prevents it from entering the Galaxy. 
Third, at the Razin frequency, the index of refraction of the 
plasma the CRs traverse suppresses the relativistic beam- 
ing of synchrotron emission. Below the Razin frequency, the 
synchrotron spectrum falls off preci pitously. In the Galaxy, 
the R azin frequency is ~ 0.1 MHz. IProtheroe fc Biermannl 
( 1996) calculated the kHz to GHz radio background, includ- 
ing synchrotron emission. In terms of energy density, the 
spectrum falls off at frequencies below about a MHz. 

The prospects for strong sub-MHz synchrotron radi- 
ation from clusters are a bit more encouraging. Galaxy 
clusters have much lower densities than galaxies, so 
bremsstrahlung and ionization losses are weaker. The low 
density and high hot temperature of the gas that fills clus- 
ters also means that free-free absorption will be much less 
effective at blocking sub-MHz radio emission (discussed fur- 
ther in § 12. 3p . Finally, the low density means that the Razin 
frequency will be much lower. For typical galaxy cluster 
densities and mag netic field strengths, the Razin cutoff is 
l|Schlickeiserll2002h 



|Pulkl|l985l). Voya q er has detected 1 
sion jKurth et all Il9841 ) 



v R = 18.5 kHz n e ^-zB^ G . 



(3) 



with n ej _ 3 = n e /(10 cm" 3 ) and B^g = B/(fiG). We may 
expect that synchrotron emission in clusters may extend all 
the way down to kHz frequencies. 

In addition, there can be low frequency synchrotron 
emission from more diffuse structures. There appear 
to be radio structures associated with poor galax; 



alaxy 
2006; 



groups or filaments of galaxies (jDelain fc Rudnickl _ 
iBrown fc R udnick 2009). Shocks from large scale structure 
formation are also predicted to produce low frequency radio 
emiss ion (|Waxman fc Loe b 2000; Kcsh etTWaxman. fc Loebl 
I2004D . 

Synchrotron emission is expected to have a brightness 
temperature limited to Tj, < 10 11-12 K, because of Inverse 
Compton scattering by the sy nchrotron emitting particles 
among other concer ns (e.g., iKellermann fc Paulinv-Tothl 
ll969l ; lReadheadlll994l 1. This would imply a low energy den- 
sity of synchrotron radio emission at low frequencies, 



I^.maxsynch < 1-3 X 10 25 ergS Cm 3 ^kHz7b,12 



(4) 



where ^ kHz = ^/kHz, and T M2 = T b /IQ 12 K. This limit only 
applies to steady-state sources; transients can have much 
greater brightness temperatures. 



2.1.2 Processes in Plasmas 

A number of plasma processes generate low frequency ra- 
dio emission. These processes are often coherent and can 
evade the Tt = 10 12 K limit on incoherent synchrotron 
emission. Many of these proc esses are observed in the Sun 
(see the review bv lDulklll985r i and the Solar System. Lang- 
muir waves in plasma can be converted into radio emis- 
sion near the plasma frequency and its first harmonic (e.g., 



4 kHz radio emis- 
n ow known to be emitted some- 
where near the heliopause (|Kurth fc Gurnettl 120031 ). Since 
its frequency is comparable to the plasma frequency of the 
ISM and the Solar Wind at the heliopause, it is gener- 
ally believed to b e powered by Langmuir wave conversio n 
(e.g., lMacek|[l996l ; iTreumann.' Macek. fc Izmodenov|["l998r i. 
Coherent plasma processes have been suggested to generate 
radio emission in AGNs, but it is n ot clear that coherent 
emission can escape from the AGNs (|Merrosejll999h . 

When the cyclotron frequency (eq. [2j exceeds the 
plasma frequency (eq. [T]), and when the electron veloc- 
ity distribution function grows with perpendicular velocity, 
the cyclotron mas er instability can also generate bright ra - 
dio emission (e.g., IWinglee fc Dulklll986l ; lTreumannll2006l ). 
The cyclotron maser mechanism is believed to be the 
main source of Jupiter's bright decam etric radio emission 
l|Hewitt. Melrose, fc Rucnnmark 1981) and Earth 's auroral 
kilometric radiation l|Wu fc Ledl 19791; IZarkaHl998T) . and may 
opera te around exoplanets (e.g., Bastian. Dulk. fc Leblanq 
21 U)i:'.. The cyclotron maser is also believed to operate in 
stars, being responsible f or the Sun's microwave spike bursts 
and stellar radio bursts (Mel rose fc Dulkll 19821 ). and coher- 
ent periodic radio emission from brown dwa rfs, though at 
higher frequencies |Hallinan et al.ll2007l . I2008T) . 

The same cosmic rays that produce the bulk of the 
synchrotron emission can also produce radiation by inter- 
acting with plasma. Electrons travelling through an inho- 
mogeneous plasma can emit transition radiation below the 
Razin frequency, cau s ed by a varying index of re f ractio n 
^Fleishman fc Kahleij Il992j ; iFleishman fc Tokarevl ll995l V 
Transition radiation is expected to be the main source of 
Gala ctic emission below ~ 100 kHz (|Fleishman fc Tokarevl 
Il995h . 



2.1.3 Exotic possibilities 

Pulsar waves - Pulsars are expected to generate magnetic 
dipole radiation of enor mous amplitude as they spin down 
l|Gunn fc Ostrikerlll969h . The pulsar wave would have a fre- 
quency equal to the pulsar spin frequency, which would be 
less than a kHz for a millisecond pulsar (MSP) and ~ 30 Hz 
for the Crab Pulsar. However, the environment around pul- 
sars is expected to be filled with dense plasma, which might 
preven t the escape of these pu lsar wave into intergalactic 
space (|Goldreich fc Julianlfl969h . 

There are r oughly 40000 MSPs in the 
Galactic disc jFerrario fc Wickramasinghel 120071 ; 
IStorv. Gonthier. fc Harding [20071; iLorimer] 120 08). and 
1000 in globular clusters IjHeinke et al. 1 120051 ). Assuming 
they have a n average lum inosity of 10 34 ergs s _1 , each last- 
ing 10 Gyr (L orimerllioOSf) , and sca ling from the stella r mass 
of the Milky Way (~5x 10 10 M ; lFlvnn et al.ll2006l ) to to- 
tal stellar mass o f t he local Universe (~ 5 x 10 8 Mq Mpc -3 ; 
ICole et all l200ll ; iHopkins fc Beacoml 120061 ). MSPs can 
provide a maximum energy density of 4 x 10 -20 ergs cm -3 
from galactic discs and 1 x 10 -21 ergs cm -3 from globular 
clusters. Young pulsars are born at approximately the 
supernova rate, or Tsiv = 1/(30 yr) for Milky Way- like 
galaxies, and have a spin energy of 5 x 10 49 erg, yielding a 
maximum ~ 30 Hz background of 6 x 10 -18 ergs cm -3 . The 
young pulsars and MSPs in the galactic discs are unlikely 
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to contribute to the extragalactic background, however, 
because of free-free absorption in the galactic discs and 
the interstellar plasma frequency is often above 1 kHz. 
Some pulsars may contribute if they are kicked out of their 
galactic disc during the supernova. Most likely, intergalactic 
absorption will reduce the energy densities further from 
these maximum estimates, so these energy densities are 
upper limits. 

Conversion of gravitational waves - Possibly one of the 
most exciting sources of low frequency radio waves would 
be gravitational waves (GWs). Typical sources of gravita- 
tional waves have frequencies of le s s tha n about a kilo- 
hertz. iMarklund. Brodin. fc Dunsbvl (2000) showed that if 
they passed through a magnetized plasma, they can gener- 
ate radio waves of the same frequency. Even if non-linear 
plasma effects generate radio waves of high er frequency, the 
resul ting frequency can still be below MHz (Marklun d et al.l 
2000). Gravitational waves can also excite cyclotron reso- 
nances l|Papadopoulosll2002l ). Black hole ringdown is pre- 
dicted to be an important source of these converted gravi- 
tational waves l|Clarkson et al.ll2004T ). 

LIGO's frequency range (~ 100 — 1000 Hz) is above 
the plasma frequency of several IGM phases, meaning that 
sources detectable by LIGO may excite radio waves. These 
sources can include compact object mergers and cosmic 
strings. Transients like merging stellar-mass black holes will 
produce bursts of gravitational waves. The cosmic energy 
density in gravitational waves is roughly u — rt hulBt u ma ^, 
where F is the rate of observed gravity wave bursts, f burst ~ 
v~ x is the duration of the burst, and « max is the gravi- 
tational wave energy density within the burst itself. This 
can be evaluated approximately as u ~ c 2 v 2 h 2 F/G, where 
h is the strain at Earth. LIGO has constrained the rate 
of h> 10" 21 Hz- 1/2 bursts with 64 Hz<^<2000 Hz to be 
Tbright < 3.6 yr" 1 l| Abbott et al.ll2009al '). implying that the 
cosmic energy density of bright gravitational wave bursts 
is w b right<4 x 10~ 16 ergs cm~ 3 ui 00 h 2 _ 2 i, where v 500 = 
^/(500 Hz) and h-21 = /t/(10" 21 Hz" 1/2 ). An associated 
radio background is likely to have a much smaller energy 
density, because only a small fraction of the gravitational 
waves is converted into radio and because intergalactic ab- 
sorption will convert the radio background into heat. How- 
ever, faint gravitational wave bursts are not as well con- 
strained by LIGO. 

LIGO has also constrained the 100 Hz stochastic grav- 
itational wave background (SGWB) to have u < 5.7 x 
10 -14 ergs cm -3 , assuming they have a po wer law spec- 
trum, which should apply to cosmic strings (jAbbott et al.l 
l2009bh . Again, any related radio background would proba- 
bly have a much lower energy density, especially considering 
intergalactic absorption. 



2.2 The Intergalactic Medium: Our Detector for 
Sub-MHz Radio Waves 

The phases of the intergalactic medium (IGM) are believed 
to fall into three main categories: a low density (<5 < 10) , 
relatively cold (~ 10 4 K) photoionized phase visible as the 
Lyman a forest; a denser, (8 ~ 5 — 100) collisionally ion- 
ized phase that is hotter (T fti 10 5 — 10 7 K), known as the 
Warm-Hot Intergalactic Medium (WHIM); and condensed 
structures (5 > 10 0) associated with galaxies and clusters 
|Dave et alJll999h . Since each of these phases have different 
temperatures and densities, they interact with radio waves 
in different ways, with different absorption properties and 
thermal contents. 

I shall consider the effects that a sub-MHz background 
would have on these IGM phases. A large enough back- 
ground would alter the properties of these phases, and would 
indirectly be in conflict with observation. I now discuss the 
phases I will be considering in further detail (properties sum- 
marised in Table I2"l2"j) : 

• The Lyman-a forest: The Lyman-a forest is detected 
as weak Lyman-a absorption lines (by definition A r (HI) < 
10 17 ' 2 cm -2 , but typically AT (HI) < 10 145 cm~ 2 ) in the spec- 
tra of bright background objects, such as quasars. The 
Lyman-a forest is particularly prominent at z>2, where 
it dominates the baryonic content of the Universe. How- 
ever, similar clouds exist at z = ; where they make up 
~ 30% of the baryonic gas (e.g., iDobrzvcki et al.l |2002| ; 
iPenton. Stocke. fc Shuul \2004 ). These clouds represent a 
'fluctuating Gunn-Peterson trough' - that is, they are denser 
regions in the background IGM, whic h by itself is too rar - 
efied to appear as absorption at low z (|Bi fc Davidsen|[l997l ) . 
Lyman-a clouds are somewhat clustered near large-scale 
structure; most do not appear to be associated with individ- 
ual galaxies, and some a re in voids (|Penton. Stocke.~fc Shull 
|2002|; IStocke et al.ll2007r i. I consider Lyman-a forest clouds 
with overdensities 8 = 0.1 and 1.0, which are typical 
of Lyman-a forest cloudjf], with temperatures of 10 4 K 
(|Bi fc Davidsenlll997h . and with radii of 100 kpc. 

• Warm-Hot Intergalactic Medium: The WHIM is one 
of the more elusive phases of the IGM, despite being pre- 
sumed to contain about w 30% of the baryoni c content 
of th e Univer se at z = (see th e review by iBregmanl 
120071 ; also see ICen fc Ostrikerlll999l ). Its presence is none 
the less inferred from UV and X-ray absorption lines of 
highly ionized s pecies, such as collisionally ionized O VI 
absorbers (e.g., iTripp. Savage, fc Jenkins 2000; Ccn et ail 
|2001| ; IDanforth et al.l 120061 ; IDanforth fc Shulll I2008T ). The 
WHIM is believed to be associated with structure formation; 
as gas gravitationally collap ses to form large-s cale structure, 
it is shocked and heats up (|Dave et al.ll200ll ). It should be 
clustered near large-scale structure but generally would not 
be associated with individual galaxies. I consider WHIM 
clouds with overdensities 8 = 5 and 100, and with a temper- 



Gravitational waves detectable by LISA, such as those 
from supermassive black hole mergers, will be below the 
plasma frequency of any reasonable density medium. They 
will be detectable only if they excite radio waves of much 
higher frequency. 



3 Note that while Lyman a clouds are 'overdense' with respect 
to the background IGM that fills most of the volume of the Uni- 
verse, they can have less than the average baryonic density of the 
Universe (8 < 1). This is because the baryonic density in most of 
the volume is much below the average baryonic density. 
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ature of 10 6 K, which are typically expected values. I assume 
WHIM clouds have radii of 1 Mpc. 

• C ool WHIM: Numerical simulations by iKang et al.l 

predict the existence of a 'cool WHIM', which is col- 
lisionally ionized and with the same densities as the WHIM, 
but colder (10 3 K < T < 10 5 K). The existence of O VI absorp- 
tion systems apparently associated with cold L yq absorbers 
may be a sign that this phase actually exists (|Tripp et al.l 
2008). I consider a cool WHIM cloud to have an overdensity 
S — 10, a temperature of 10 4 K, and a radius of 1 Mpc. 

• Cluster gas: An extremely hot (10 8 K) phase of gas 
pervades galaxy clusters, and is easily visible in free-free 
X-ray emission. I consider galaxy clusters with gas den- 
sities of 10" 4 cm -3 (2 x 10 12 M Mpc -3 or 5 = 400) 

-3 



and 10~ 3 cm~ 3 (2 x 10 13 M M pc" 3 or 8 = 4000) 
|Mohr. Mathiesen. fc Evrardl 119991 ). The latter density is 
more typical of rich clusters. Cluster g as temperature is as- 
sumed to be 10 8 K (D avid et al.lll99a l. and the cluster gas 
has a radius of 1 Mpc in this paper. 

• Metal-line absorbers: Relatively dense photoionized gas 
is also present in the inner tens of kiloparsecs around L* 
galaxies, manifesting itself as strong Lyct and metal ab- 
sorption lines in quasar spectra. I consider strong Mg II 
absorbers wi th n = 0.01 cm" 3 and R = 1 k pc with 
T = 10 4 K (|Ding et all 120031 ; lEllison et all I2004T ). Weak 
metal lin e absorbers are also present around sma ller, dwarf 
galaxies |Stocke et al.ll2004l ; IKeenev et al.| [2006). and pos- 
sibly in intergalactic space. Photoionized weak metal-line 
absorbers include dense structures with n « 0.01 cm -3 
and R ~ 10 pc, detected in low ionization lines like Mg II 
and F e II (|Rigbv. Charlton, fc Churchillll2002l ; iRichter et all 
2009); kiloparsec-scale structures with n ~ 10 cm that 

are observed in C IV at z ~ 1 and Mg II at z « 

llRauch. Sargent, fc Barlowl l200ll; ICharlton etHI 120031 ; 



iMilutinovic et a" 20061 ; Narayanan et al. 20051 ); and possi- 
bly high-ionization absorbers with n ~ 10~ 5 cm~ 3 and R ~ 
100 kpc seen in high ionization lines, such as C IV haloes 
at z k, (|Chen. Lanzetta. fc Webbll200ll ; lLacki fc Charlton 
2010) and intergalactic O VI clouds at z ~ 1 (IZonak et al 



2004; lTripp et al.ll2008l ; lOppenheimer fc Davel l2009). I con- 
sider examples of all these structures, with a typical pho- 
toionized temperature of T — W 4 K. 



an IGM cloud-detector would be expected to last. There are 
four relevant time-scales for most IGM phasefl: 

(i) Dynamical time-scale - The time for the cloud to 
collapse under self-gravity. I assume a dynamical time of 
idyn ~ l/v 7 ^, or 

idyn « 190 Gyr S~ 1/2 (5) 

(ii) Thermal time-scale - The IGM can cool radiatively, 
but is also being photoionized by the extragalactic UV back- 
ground. We can define a thermal time-scale 

3kT 



^therm 



(6) 



In the IGM, the cooling coefficient Ajv can be negative be- 
cause of the ph otoionizing background. I have use d the val- 
ues of Ajv from I Wiersma. Schave. fc Smith! (|2009l ). 

(iii) Sound crossing time-scale - In general, this is not 
equal to the the dynamical time-scale, since IGM clouds like 
those in the Ly man-a forest are not virialized. The sound 
speed is c s — */ 10fcT/ (3m.tr), giving a typical sound cross- 
ing time of 

1/2 



59 Gyr R Mp cT~ 



(7) 



where Rm P c is the radius of the cloud in Mpc, and T4 = 
T/10 4 K. 

(iv) Age of the Universe - For simplicity, I use tu ~ 
10 Gyr. 

I choose the shortest of these times to be £igm , the expected 
survival time of an IGM cloud. These times are compared 
in Table [2~2l 

Finally, there is the relativistic phase of the IGM, the 
cosmic rays (CRs) . Cosmic ray electrons (and positrons) are 
known to exist in galaxy clusters, through their synchrotron 
radio emission and the Inverse Compton emission in X-rays 
l|Ferrari et al.ll2008l ; [Re"phaeli et al.ll2008l ). Presumably there 
also is a more diffuse cosmic ray background that pervades 
the Universe. A large fraction of cosmic rays, both elec- 
trons (and posi trons) and protons, escap e from smaller, less- 
dense galaxies (|Aublin fc Parizoti l2006). Dark matter anni- 
hilation or decay may also inject CRs into the IGM (see 
IChen fc Kamionkowskii r2004). However, no truly diffuse ex- 
tragalactic CR background has ever been detected, except 
at the highest energies. 



These phases span a wide range of density, tempera- 
ture, and locations near structures; from Lyman a clouds 
that often appear in voids, to metal line absorbers that are 
concentrated very near galaxies. However, most of the vol- 
ume of the Universe is filled with voids with extremely low 
density gas (5 <C 1). Thus, if there was sub-MHz radio that 
somehow filled only voids, then it may not have a great ef- 
fect on detected IGM phases. Even though there are some 
Lyman-a clouds in voids, they would have a plasma fre- 
quency greater than the background IGM; thus they would 
reflect incident radio waves of the lowest frequencies rather 
than being directly heated by them. Even so, the pressure 
exerted by the exterior radiation background on the cloud 
could still have an effect if it was too high. 

The sensitivity reached by an IGM cloud-detector de- 
pends on its effective integration time, or the time tiGM that 



2.3 Propagation of Sub-MHz Radio through the 
IGM 

Sub-MHz radio waves suffer a number of absorption and 
scattering processes as they traverse the IGM. 

Thomson scattering - Since the IGM is ionized, it con- 
tains free electrons that will Thomson scatter incident elec- 
tromagnetic radiation of all radio frequencies. However, the 
optical depth to Thomson scattering is very low, since the 
IGM is so rarefied: 

4 There are additional possible time scales in some phases. For 
cluster gas, there is also the Alfven crossing time, tAifven = R/ua 
where da = B/^A-rrp is the Alfven speed. The Alfven crossing 

time is tAifven = 14 GyriiMpcB^Q^^ » w hich is longer than the 
sound crossing time for clusters. 
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Table 1. IGM Phase Summary. 



Phase 


R 


n e [S e ] 


Vp 


T 


"at 


^•therm 


^dyn 


^sound 


*IGM 




(Mpc) 


(cm -3 ) 


(Hz) 


(K) 


(erg cm 3 ) 


(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 


Background lOlVl 


1000 5.0 x lO" 10 [0.002] 


0.2 


i n4 
10 






4000 


6 x 10 


10 


Underdense Lya forest 


0.1 


2.5 x 10" 8 [0.1] 


1.4 


10 


— 2.3 X 10 


200 


600 


6 


6 


Lya forest 


0.1 


2.5 x 10~ 7 [1] 


4.5 


10 4 


-1.7 x 10~ 23 


30 


200 


6 


6 


Low density WHIM 


1 


1.25 x 10- 6 [5] 


10 


~ 10 6 


1.4 x 10~ 23 


800 


80 


6 


6 


High density WHIM 


1 


2.5 x 10~ 5 [100] 


15 


~ 10 6 


2.1 x 10~ 23 


30 


20 


6 


6 


Cool WHIM 6 


1 


10- 6 [4] 


9.0 


10 4 


-1.5 x 10~ 23 


9 


90 


60 


9 


Low density cluster 


1 


10~ 4 [400] 


90 


10 8 


3.7 x lO" 23 


400 


9 


0.6 


0.6 


High density cluster 


1 


10- 3 [4000] 


280 


10 s 


3.0 x 10~ 23 


40 


3 


0.6 


0.6 


Strong Mg II absorber 


0.001 


10~ 2 [40000] 


900 


If) 4 


-6.5 x 10~ 25 


0.02 


1 


0.06 


0.02 


Dense weak Mg II absorber 


10" 5 


10" 2 [40000] 


900 


10 4 


-6.5 x 10~ 25 


0.02 


1 


6 x 10" 4 


6 x 10~ 4 


Weak Mg II/C IV absorber 


0.001 


10~ 3 [4000] 


280 


10 4 


-1.2 x 10~ 24 


0.07 


3 


0.06 


0.06 


C IV/O VI Galaxy Halo 6 


0.1 


10~ 5 [40] 


28 


10 4 


-1.1 x 10~ 23 


1 


30 


6 


1 


a Cooling rates for intergalactic clouds at 


z = for Z = 


3.1Zrr\ metallicitv fromlWiersma 


et al.l ((2009 


). Neg 


ative cooling 


rates 



imply net heating. 
These phases have not been confirmed observationally. 




Figure 1. Absorption path lengths for the Coma Cluster at low 



frequencies, with B = 0.5 £tG and n 



10" 



The thick 



black line is free-free absorption. The thin diagonal lines are for 
synchrotron absorption; solid assumes the spectrum falls precipi- 
tously at low energies (yc ^ 30 MHz) , long dashed assumes a flat 
spectrum at low energies, short dashed assumes an E~ 2 spec- 
trum continuing down to 7 = 1, and dotted assumes an £? -3 
spectrum continuing down to 7 = 1. The horizontal line is Thom- 
son scattering. Synchrotron absorption dominates at frequencies 
above the Razin cutoff (vr), but free-free absorption remains be- 
low that frequency. The cluster is optically thick at v < kHz from 
free-free absorption, and near the Razin cutoff from synchrotron 
absorption if the low energy CR electron spectrum is steep. 



1950 5' 1 Gpc. 



(8) 



Therefore, Thomson scattering is unlikely to be of any sig- 
nificance in the IGM. 

Free-free absorption - Free-free absorption becomes ef- 
fective at lower frequencies, and would blanket out the low 
frequency Universe. The free-free absorption coefficient is 



0.018T" 3/2 giFn^" 2 cm 5 Hz 2 K 3/2 . 



(9) 



The Gaunt factor gg for radio frequencies above the plasm a 
frequency and below kT/h is (|Scheueijll960l ; ISpitzedfl97Sl ) 

gs = 9.77[1 + 0.1951og 10 T 4 - 0.1301og 10 i/M Hz ]. (10) 

Over the frequencies considered in this paper, ~g~ft ~ 10 — 20. 

The absorption mean free path of sub-MHz radio in the 
IGM is roughly 



A a = 19 Mpc T! /2 gWJE 



-2 2 

^kHz 



(11) 



The background IGM phase I consider (<5 e = 0.002) is 
transparent across Gpc distances down to ~ 15 Hz. How- 
ever, a radio wave could run into denser structures that are 
more opaque and be absorbed. Typical Lyman-a and WHIM 
clouds become transparent above about 10 — 1000 Hz, as- 
suming scattering within the clouds is negligible. The strong 
Mg II absorbers are opaque for v < 250 kHz, even without 
scattering. Thus, the time \ a /c that a radio photon can 
travel is limited by the time until it reaches one of these 
denser structures. Somewhat perversely, scattering can ac- 
tually increase the absorption time, by trapping photons in 
IGM phases with low absorption. 

Synchrotron absorption - Synchrotron absorption 
would be effective in galaxy clusters, where a high content 
of CRs and magnetic fields are present. The synchrotron 
absorption coefficient is 



8tt^ 2 



dEP s 



N(E) 



£ 2 



(12) 



where P s (y, E) is the synchrotron power at frequency u 
(above the Razin cutoff; equation [3} of a single elec- 
tron/positron at energy E: 



P s (y,E) = 

with 



sin ct 



3j 2 eB sin a 
47rmc 



(13) 



(14) 



The synchrotron absorption depends on the CR elec- 
tron spectrum in clusters. It is possible to estimate the elec- 
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tron spectrum for 10 3 < 7 < 10 4 from the MHz to GHz syn- 
chrotron radio they emit using the monochromatic approx- 
imation (|Felten fc Morrison! Il966l ; ISchlickeiserll2002l ): 



„ esma 1 

N(E) = 9— i -57 £„. 

m 2 a T c 4 B 



(15) 



The disadvantage of this method is that there is a minimum 
CR energy probed by synchrotron radio emission, because 
there is a minimum radio frequency observed v m i no b s : 



B \ ~ 1//2 



1/2 



(16) 



We can, of course, use models for the spectrum at lower 
energy, but the presence of sub-MHz radio would actually 
alter the low energy spectrum through synchrotron heating 
(as described in § [BJ. When possible, we should use clus- 
ters with fluxes detected at the lowest possible frequencies. 
The Coma cluster radio halo has observations at frequencies 
ran ging from 30 MHz to 5 GHz, so I have us ed it as an exam- 
ple (|Thierbach, Klein, fc Wielebinskil \200$) . Other clusters 
are typically observed at GHz frequencies, though more clus- 
ters will be seen at lower frequencies (down to ~ 15 MHz) 
with LOFAR. 

In Figure [T] I show the importance of Thomson scatter- 
ing (grey), free-free absorption (red), and synchrotron ab- 
sorption (blue) for the Coma Cluster at low frequencies. 
The synchrotron absorption depends sensitively on the low 
energy (7 < 1000) CR electron/positron spectrum. I show 
that if there are no electrons below 7minobs (solid), syn- 
chrotron absorption is relatively weak at lower frequencies. 
The steeper the spectrum below 7 m inobs, the greater the syn- 
chrotron absorption, as can be seen if the spectrum has a 
constant dN/dE for 7 < 7 m i no bs (dotted) and if the spectrum 
goes as dN/dE oc E~ 2 or E~ 3 for 7 < 7 m i n oba (dashed). 

Plasma inhomogeneity scattering - Finally, sub-MHz 
radio may scatter off pl asma inhomogeneities in the IGM 
(|Cohen fc Cronvnl [l974). Suppose the plasma inhomo- 
geneities have a spectrum P n (q) = C^q~ a , where the 
wavenumber q — 2n/£ ranges from a minimum go (large 
wavelength) to a maximum qi (small wavelength). The 
power P n is defined by J P n (q)d 3 q — (An 2 ). For an isotropic 
Kolmogorov spectrum of inhomogeneities (a = 11/3), where 
£0 3> £\, the coefficient C 2 is 



C' n 



,1/3 



(An 2 ) 



37r 4 /3 £2/3 



(17) 



The image of a point source is smeared over an angle 
= c/ (2ttv£ c ), where 



An 2 -^sC 2 f(a) 



-V(a-a) 



where s is the path length through the medium, and 
l2 _ a T(a/2)(a-2) 



/(<*) = r 



r(2-(o/2)) 



(18) 



(19) 



but equation [181 holds only if 2 < a < 4 and £\ < 4 < l o 
l|Lee fc JokipiTll975l ; ICordes. Weisberg. fc Boriakoflj Il985l ). 
If we take the scattering mean free path A s to be the dis- 
tance when the image of a point source is smeared out over 
8 C = 1, then £ c = c/(n 2 u). For a Kolmogorov spectrum of 



turbulence (a = 11/3), 



A 



i 1.5 x 10- 5 Mpc ^Mpc^n" Vm£ (2°) 
« 240 Mpc e 2 % pc A^67 2 ^Z (21) 

where A n = (An 2 )/n 2 . The scattering coefficient is then 
just defined to be <5„ = A^ 1 - 

There is a great deal of uncertainty in the derived 
path length. Most importantly, equation [20] holds only if 
£1 < cI(-k 2 v) « lo- I n the Milky Way, h « 10 7 cm 
jSpangler fc Gwinnl [l99fj| ) ; if this is true in the IGM, then 
these formulae are valid only if v <C 300 Hz. The two natu- 
ral scales fo r £1 are the ion Larmor r adius and the ion iner- 
tial length (|Spangler fc Gwinnlll990l ). both of which are ex- 
pected to be very large in the IGM. The mean free path also 
depends on the outer scale of the largest inhomogeneities 
£0 and their amplitude (An 2 ). In our Galaxy, the inhomo- 
geneities continue m erely up to ~ pc scales rather than the 
size of the Galaxy (|Haverkorn et al1 [2008). On the other 
hand, the amplitude of the fluctuations is almost certainly 
less than the mean ISM density in our Galaxy. 

The presence of scattering in the IGM would trap sub- 
MHz radio near their sources of emission. The scattering 
optical depth of a medium with radius R — i?M pc Mpc is 
t s = R/\ s , or, if the estimate in equation 1201 holds. 



T s W 6.5 X 10"7?Mpcl Oi Mpc A ™ n e,-3^MH: 



11/3 



4.1 x 10 s i4 



-2/3 



A„<5 e i/ kHz 



11/3 



(22) 
(23) 



tM P c *0,Mpc 

In an optically thick medium, the sub-MHz radio must dif- 
fuse with a typical time scale idiff = 3-R 2 /(cA s ), or 



tdiS 



6.4 x 10 11 yr R 2 



2/3 

Mpc^o.Mpc^- 1 " 



A 2 -11/3 



4.0 x 10 15 yr R 2 



-2/3 



Mpc'o^pc^"^ "kHz 



11/3 



(24) 
(25) 

Thus, based on the short distance found in equation 1201 
and long diffusion times, even if we could somehow place a 
sub-MHz detector outside of the Galactic WIM, we still may 
not see most of the sub-MHz radio sources in the Universe. 
However, since this is not an issue at present, the scattering 
can serve to our advantage. Sub-MHz radio is trapped near 
its source; therefore, if we observe the effects of sub-MHz 
radio (as described in following sections), we know that its 
source is nearby. Trapping localises sub-MHz radio for us, 
and we do not simply have to consider a diffuse background. 

It is possible to define an effective absorption mean free 
path, which is the displacement a radio photon can travel be- 
fore being absorbed. When scattering is strong, the effective 
absorption mean free path is much shorter than X a , because 
the scattering traps the photon near its source. The effec- 
tive absorption mean free path is A aj0 ff = 1/ yj ct v (ct v + 8 V ). 
If equation [20] holds and free-free absorption is the main 
absorption source, then ^ > q„, and 



l/2 T 3/4 17/6 



(26) 



A a , cff = 6.9X10 4 Mpc£l% pc A~ 1/2 S: 2 gWJE~ 1/2 T! 
We see that A aiC ft steeply falls as the frequency decreases. 



3 LIMIT FROM Q R 

Since the extragalactic ultra-low frequency radio back- 
ground cannot be detected from the Solar System, it can 
be thought of as dark radiation, a relativistic component of 
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the Universe not directly seen. Like other postulated kinds 
of dark radiation, sub-MHz radio must have a gravitational 
effect on the expansion of the Universe. So cosmological 
bounds on the radiation energy density of the Universe also 
limit extragalactic sub-MHz radio. 

Most cosmological limits are very stringent for the 
early Universe (z > 1000), but do not rule out the growth 
of Qr late in the Universe. Since presumably sub-MHz 
radio would be radiated by non-thermal processes after 
the CMB has decoupled and structure has begun form- 
ing, limits from the CMB and Big Bang Nucleosynthesis 
are largely irrelevant. Type la supernovae data give us the 
limit Q r<0.2 for z < 1000, or it < 1.7 x 10" 9 erg cm" 3 at 
z = (Ze ntner fc Walkerl 1200 j ). This argument applies to 
the globally averaged sub-MHz radio background, and not 
to specific, compact regions like clusters. 

Unlike other more exotic forms of dark radiation, sub- 
MHz radio does interact with intergalactic matter. These 
interactions are then visible from Earth, and are generally 
much more constraining. 



4 IGM HEATING FROM FREE-FREE 
ABSORPTION 

Just as free-free absorption prevents sub-MHz radio from 
entering the Galaxy, it should also be present in intergalac- 
tic space. The free-free absorption converts low frequency 
radio waves into IGM h eat, which can affect the dynamical 
ilCiotti fc Ostrikerj |2004 ) or thermal state of the gas. There- 
fore, the IGM serves as a low-frequency radio detector: if 
the radio background was too high, the IGM would heat up 
to a higher temperature than observed. Low frequency ra- 
dio emission would create anomalous heating in the IGM; 
the lack of such heating puts an upper limit on the radio 
background. 



4.1 Energy density limits 

Consider a small sphere of IGM plasma immersed in an 
isotropic bath of sub-MHz radio with steady intensity J v for 
a time £igm (§ 12. 2[) . during which we expect the IGM not 
to be otherwise heated or cooled significantly. The heating 
rate of the plasma parcel is 



Q w ttA£ 



-d.J v 
d£ 



dv. 



(27) 



where A is the surface area of the sphere, t is the mean path 
length through the sphere, and ^jf- is the mean change in 
radio intensity per path length. If the brightness tempera- 
ture of the radio bath is much greater than the temperature 
of the plasma, then =4jf- w — «„ J„, given an absorption coef- 
ficient ol v . We shall suppose that we are considering a very 
small parcel, so that there is no scattering and I ~ R, the 
radius of the sphere; then At w 3V , where V is the volume 
of the sphere. To avoid excess heating, 



Q < 



SNkTiGM 

2tlGM 



(28) 



In this equation, N is the total number of particles and is 
approximately equal to twice N e , the number of electrons. 



Equation [27] applies to a small parcel of gas, and the 
relevant J„ is the intensity within that parcel. However, the 
intensity of the radio bath may vary drastically throughout 
the cloud if it is externally illuminated and optically thick 
to free-free absorption, in which case the radio intensity is 
very low deep inside the cloud. 

Spatially homogeneous J„ case - First, consider the case 
when the radio background is constant throughout an IGM 
cloud. This assumption is appropriate if either the cloud 
is optically thin (to both scattering and absorption), or if 
internal objects distributed nearly continuoushjj within the 
cloud are the source of the radio bath, as would be the case 
for synchrotron or transition emission emitted by CRs in the 
cloud. 

Suppose the radio bath consists of a nearly constant 
bump at frequency v and with width Av<v. Then condi- 
tions [27] and [28] gives us 



27rflGM 



(29) 



Assuming that n e — 8 e {nb) = n/2, the free-free absorption 
coefficient gives the following upper bounds: 



Ai/J„ < 2.1 x 10 



T 5/2 2 

Avu v < 8.7 x 10 



T 5/2 2 
-11 1 4 y kEz 

ilO<7ff,15<Se 

2 

kHz 



-2 -1 

ergs cm s sr 



iio<7ff,i5<5c 



erg cm 



1 (30) 
(31) 



with T 4 = T/(10 4 K), */ kHz = z//(kHz), gWJE = ffff/15, and 
tw = t IG M/(10 Gyr). 

In Table [2] I consider various phases of the IGM and set 
limits on the homogeneous radio background within these 
phases. The best limits come from dense, cold phases of the 
IGM, such as galactic metal-line absorbers. Hot phases like 
the WHIM and cluster gas place relatively poor limits on 
the radio background. The corresponding brightness tem- 
peratures of the allowed radio backgrounds are high enough 
that ignoring the thermal emission of the IGM cloud is ap- 
propriate. 

Although I assumed the spectral intensity was a bump 
of width A^< v when calculating the limits in equations 1301 
and 13 II the heating can be integrated over any desired spec- 
trum in equation 1271 In particular, for a power law of in- 
tensity J„ = Jo{v/vo) q in the frequency range i/ m m <C v <C 



Jo < 2.1x10" 



n5/2 



Setwga, 15(^0 



^max 
9-0-1 

(l.l-9)?Nr 



vq In 



(ff>l.l) 
(«<1.1) 
(9 = 1-1), 



where Jo is in cgs units (ergs cm -2 s _1 sr -1 ) and I have 
used the approximation that the Gaunt factor goes as v~ 0A 
with a normalisation gs z iz(yo) = flff(wo)/15 at v$. 

Externally illuminated, optically thick case - Now sup- 
pose an optically thick cloud has no internal sub-MHz lu- 
minosity, but is immersed in an external radio bath. Since 
the cloud is optically thick, the interior of the cloud will at 

5 Continuous means here that the separation between sources is 
less than the effective absorption mean free path within the cloud. 
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first be unaffected by the radio bath. The exterior of the 
cloud will heat up through free- free absorption. As it heats 
up, it will become more transparent to free-free absorption 
(equation [9]), both because of its higher temperature and 
because it can expand away and dilute. The next layer in is 
then exposed to the external radio bath, so that the cloud 
is effectively eroded by the external radio bath. 

I assume that the intensity of the radio bath is homoge- 
neous and equal to its external value J„ within one absorp- 
tion effective mean free path A aje fl = [ct v ,a{oii/,S + 8 V )]~ 1 / 2 
of the cloud exterior, where 8 U = A^ 1 is the scattering 
coefficient. Then the time for this layer to evaporate is 



kevap (A a ,eff 



£evap(Aa 5 eff ) 



Q/Q, or from equation [ 
3n e fcT 



ira v> ff(AvJ u ) 



(33) 



assuming the spectrum of the radio bath is a narrow bump 
of width Av. 

The time for the entire cloud to evaporate is therefore 

tevap(i?) ~ (i?/A aj eff )t e vap(A a ,eff ), Or 



>(*) 



3Rn e kT 



'1+-^-. 



(34) 



Note that there is now no frequency dependence if there is 
no scattering (<5„ = 0). Since the cloud is optically thick, it 
absorbs all of the incident radiation, which is converted into 
heat. Scattering slows down the evaporation process, since 
the incident radiation cannot penetrate as deeply into the 
cloud at any given time. 

In order for the cloud to survive in the radio bath, we 
require that t cvap (-R) > tiGM- For typical values in the IGM, 
the limits on external radio backgrounds for optically thick 
clouds ar«Q 



AT <T on i n-12 5 e TaRwpc L , 8 V 

AvJ v < 3.2x10 — 4 /In- 



tro 



OLv.fi 



Avu v < lAxW- 2l SeTiRMpc Jl + ^, 

tlO V a v> s 



(35) 
(36) 



where AvJ v and Avu v are in cgs units (ergs cm"" s~ 
and ergs cm -3 , respectively). 

The optical depth of a cloud is greater at lower fre- 
quency. Therefore, the homogeneous, optically thin case 
(eg. 1301 and eqn:FFBumpLimitu) is appropriate for the high 
frequency radio background, while the optically thick case 
(eq. [35] and I36|) is appropriate for the external lowest fre- 
quency radio background. When there is no scattering, the 
transition between the two cases occurs when rg = a v ^R = 
1; for typical values in the IGM, this is 



VcAR) « 230 Hz Z 



3/2 



£ 1/2 D 1 / 2 



(37) 



6 Note that the bounds on AvJ v and C\vu v in the externally 
illuminated, optically thick case (eq. 1351 and I36H do not match 
the homogeneous case (eq. I30l and l31t when r a c g = 1 (at either 
f ca l(i?) given in eq. 1371 or at v c3 x S (R) given in eq. 1 38 P - The dis- 
crepancy is caused by a geometrical factor: in the homogeneous 
case, the heated volume is a sphere with radius R, but the op- 
tically thick case assumes the heated volume is a thin shell with 
radius R and depth A a e g , which will overestimate the volume if 
A a , e ff = R- 



The homogeneous approximation should be valid down to 
below a kHz for typical IGM clouds, if there is no scatter- 
ing. With scattering, the cloud becomes effectively optically 
thick when r e ff « U a v ^8 v R = 1. If we naively apply the ef- 
fective mean free path from equation l26l we get the following 
characteristic frequency 



fcai.sOR) « 20 kHz T 4 



-9/34 cl2/17. 



ffff,15 



3/17 A 3/17p6/17„-2/17 



''Mpc 0,Mpc 

where A n = {An 2 )/n 2 . Below this characteristic frequency, 
the upper bound on AvJ v would go as -J 8 v ja v oc v~ h ^ . 
However, as explained in § 12.31 the approximations used in 
5v are unlikely to be appropriate for large scattering optical 
depths. 

However, I have ignored the time it takes for radio waves 
to diffuse into the centre of the cloud when scattering is 
present. Comparison with the diffusion times derived in § 12.3 
shows that the diffusion times are likely to matter at low fre- 
quencies if this scattering is present. In the case with scatter- 
ing, we must require that the radio waves have enough time 
to diffuse through the cloud in order to heat it. Conserva- 
tively, I will choose the requirement that the time to diffuse 
through the cloud radius R is less than or equal to £igm- 
From equation 1241 this gives us a characteristic frequency: 



(38) 



fdifl 



34 kHz t 



-3/11 p 6/ll 



10 



-2/11 A 3/llc6/ll 



(39) 

where tio = iiGM/(10 Gyr). If scattering is present, the free- 
free absorption bounds do not apply below this frequency. 
There may be a way around this problem - the outer layers 
of the cloud may disperse as they 'evaporate', allowing new 
radio waves to free stream into the deeper layers of the cloud. 

In Table [2j I have calculated these limits for the various 
IGM phases considered here. I have assumed A n = 1 and 
£0 — R for the case with scattering. As before, the Lyman-a 
forest limits are the strongest, because of the low thermal 
energy density of the Lyman-a forest and its long lifetimes. 
Furthermore, the low density of the Lyman-a forest means 
that it is more transparent, allowing more of the cloud inte- 
rior to be heated. Scattering weakens the limits considerably 
for i/< fcal,s(-R), since the cloud becomes effectively opaque 
at much higher frequencies than without scattering, shield- 
ing the interior from the external radio bath. Requiring that 
the diffusion time through the cloud be less than tiGM weak- 
ens the limits even more, since then there are no free-free 
absorption bounds below i><ns. For 6 = 1 Lyman-a forest 
clouds, this allows large external radio backgrounds below 
13 kHz. 

The estimates in this section ignore another consider- 
ation that may weaken these bounds. I have assumed that 
once a layer of the cloud is heated, it disperses and no longer 
contributes to the opacity. While their opacity is expected 
to go down as they heat up, the outer 'evaporated' layers of 
the cloud may still be sufficiently opaque to shield the inte- 
rior of the cloud from the external radio background. This 
will slow down further evaporation when t> 1, which will 
especially be a problem when scattering is strong. 

In this section I have only considered bounds from the 
heating of the entire cloud. There may other, more subtle 
statements that can be made. For example, even if the sub- 
MHz radio cannot penetrate deep into the interior of the 
cloud, it can still heat the exterior of the cloud, and cause 
it to become more highly ionized. The heated layers would 
then be visible in absorption lines like C IV or O VI at 
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Table 2. IGM Free- Free Absorption Upper Limits on the sub-MHz Background 



Phase 


Vp 




fcal,s(R) i> 


dm 




Upper Bounds on AfM„ 












Homogeneous (MHz) d 


Opaque (No scattering) 6 


Opaque (Scattering)-^ 




(Hz) 


(Hz) 


(kHz) 


(kHz) 


(erg cm" 3 ) 


(erg cm -3 ) 


(erg cm" 3 ) 


Background IGM 


0.20 


15 


1.2 


14 


6.7 X 10 


2.8 X 10 


2.0 X 10 9 


Underdense Lya forest 


1. 1 


7.5 


2.2 


1.8 


2.2 X 10 -13 


2.4 X 10 -23 


2.2 x 10 -18 9 


Lya forest 


4.5 


72 


11 


17 


2.2 X 10" 14 


2.4 x 10" 22 


5.5 x 10" 18 9 


Low density WHIM 


10 


41 


18 


93 


3.2 x 10~ 10 


1.2 x 10" 18 


3.0 x 10" 13 9 


High density WHIM 


15 


780 


150 


480 


1.6 x lO" 11 


2.4 x 10~ 17 


1.0 x 10" 12 9 


Cool WHIM 


9.0 


870 


50 


74 


3.7 x 10" 15 


6.2 x 10" 21 


1.9 x 10" 17 9 


Low density cluster 


90 


110 


120 


1900 


3.1 x 10" 6 


9.5 x 10" 14 


1.7 x 10" 7 9 


High density cluster 


280 


1100 


620 


6700 


3.1 x 10~ 7 


9.5 x 10~ 13 


4.3 x 10" 7 9 


Strong Mg II absorber 


900 


2.4 x 10 5 


6200 


4800 


1.7 x 10~ 16 


2.7 x 10" 17 


1.3 x 10~ 14 


Dense weak Mg II absorber 


900 


26000 


2100 


2400 


5.5 x 10" 15 


9.5 x 10" 18 


5.1 x 10" 14 


Weak Mg II/C IV absorber 


280 


26000 


1300 


1000 


5.5 x 10" 16 


9.5 x 10 -19 


1.8 x 10" 15 


C IV/O VI Galaxy Halo 


28 


2700 


150 


200 


3.3 x 10~ 15 


4.9 x 10" 20 


1.4 x 10" 16 



a Frequency at which the cloud becomes optically thick to free-free absorption, assuming no scattering. 

b Frequency at which the cloud becomes effectively optically thick to free-free absorption, assuming the scattering mean free path 
is given by eg. 1201 Note the kHz units. Assumes £o = R and A n = 1. 

c Frequency at which the diffusion time through the cloud equals iiGMi assuming the scattering mean free path is given by eg. 1201 
Note the kHz units. Assumes £q = R and A n = 1. 

d Upper bound at MHz if the radio is distributed evenly through the cloud. The upper bound on Auu v roughly scales as u 2 , with 
some logarithmic dependence from gTf . 

e Upper bound on incident external radiation with v < v ca .i(R) from cloud evaporation, if the clouds do not scatter radio waves. 
f Upper bound on incident external radiation at v = v ca \ S (R) from cloud evaporation, if the cloud's scattering mean free path is 
given by eo. 1201 and if the radio has enough time to diffuse through the entire cloud (applies if i^dift < ^cal.s(-R)). Below ^ C al,s(-R)j 
the energy density bound is weaker, going as t/~ 5 / e . Eg. 1201 may not be valid at large optical depths. Assumes £q = R and A n = 1. 
9 Since i^diff > v CB \ s(-R)i the radio may not have enough time to diffuse through the entire cloud to evaporate it when u = f C al,s(R)- 
Therefore these limits may not apply. 



the same velocity as the cloud interior. The number statis- 
tics of highly ionized absorbers compared to lower ionization 
absorbers at the same velocity may then set limits on the 
sub-MHz background even in the case of strong scattering 
and high optical depth. 



4.2 Luminosity Density Limits 

We can also set limits on the luminosity density of low fre- 
quency radio sources within the IGM. The energy density of 
a region is simply the luminosity in that region integrated 
over the time the photons could have lasted in the IGM: 

His — (Lutv . 

Homogeneous case - In the case when the sources are 
distributed inside the cloud continuously (equation I29|) . the 
limit on the luminosity density is 



2nkT 



(40) 



Suppose that free-free absorption is the main source of ab- 
sorption, with synchrotron absorption negligible. As dis- 
cussed in § 12.31 this is probably a good assumption for 
most of the IGM; even in clusters, it will hold for radia- 
tion with frequencies below the Razin frequency. If we ig- 
nore other sources of absorption, ct v = ct Vl g. We can define 
an absorption time-scale for that particular parcel of gas as 
tabs = (ca„) -1 , which is the time-scale that a photon of 
frequency v that is trapped within that gas parcel will be 
absorbed. This gives us 



2nkT tabs 

ilGM t u 



(41) 



The choice of t v is potentially very complicated, since 
it is spatially averaged. If the photons are trapped in a 
homogeneous IGM phase, t v = min[(t a bs), £igm], where 
(tabs) = {caj -1 is the mean absorption time-scale for a 
photon of frequency v over the path it travelled. 

In the spatially homogeneous opaque limit with only 
free-free absorption, when t v = (t a bs) = tabs = tff, we enter 
a simple 'calorimeter' limit: 



2nkTiGM 



tlGM 

which can be evaluated a^l 
Ave v <3A x lO 5 L Mpc" 3 (5 e r4 



tlGM 

Gy7 



(42) 



(43) 



This equation simply says that all of the radio emission in 
a region heats up the local IGM, and the luminosity inte- 
grated over time cannot exceed the thermal content of the 
IGM. The ideal properties of an IGM parcel as a radio de- 
tector are obvious: it must be relatively low density and 
cold, and it must last a long time. However, the frequency 
at which the IGM becomes a radio calorimeter increases 
with density. Low density IGM phases are good detectors at 
low frequency, while high density IGM phases are better for 
high frequency. In the case when other sources of absorption 
are present, the true limit of Ai/e v is larger by a factor of 
iff /tabs, since there are other sinks for the luminosity. 

The exact frequency f ca i at which an IGM cloud be- 
comes calorimetric depends on the radiative transfer within 



Note that 1 erg cm 



7.7 x lO 39 L Mpc" 
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the cloud. The longest possible time a radio photon can have 
to be absorbed (thus heating the cloud) is £igm, since the 
cloud has a finite age. This case applies when scattering 
traps radiation effectively. The absorption time equals *igm 
at the frequency 



^cai(ctiGM) = 13 kHz S e gs^ 5 1/2 tY 2 T 4 3/4 



(44) 



The shortest time a radio photon can have to cross the cloud 
and escape is R/c, which applies when radio emission free- 
streams out of the cloud. The frequency when the absorption 
time is R/c is simply the frequency at which it would become 
opaque with no scattering, f ca i(P), given in equation 1371 I 
calculate the frequencies when ts equals each of these times 
in Table El 

The calorimetric limit applies for all frequencies vp < 
v < VcaX, where u cal is between ^(tiGM) and v cai [R/c). In 
Table [3] I calculate the calorimetric limits on the luminosity 
density of sub-MHz radio for each phase, assuming free-free 
absorption is the only source of absorption. The Lyman- 
q forest sets very good limits on the luminosity density of 
the sub-MHz Universe. The existence of 8 = 1 Lyman-Q 
clouds implies that Aue v <6 x 10 4 Lq Mpc -3 for frequen- 
cies Vp < v < i/ ca i. For comparison, this is smaller than the 
gamm a-ray emissivity of the Universe (|Coppi k, Aharonianl 
1 19971 ). and only a tiny fraction of the bolometric emissiv- 
ity of the Universe. Denser and hotter clouds, like WHIM 
clouds or metal line absorbers set relatively poor limits on 
the sub-MHz luminosity density. 

Optically thin case - If the Universe is optically thin 
to a radio photon, so that the mean absorption time-scale 
along its path length is much longer than the age of the 
Universe, however, then t v = tn since the Universe has been 
filling up with radio emission for its entire history. Note that 
the condition (t a b s ) ^ tu must account for the fact that 
the radio photon may traverse different IGM phases with 
different local absorption time-scales i a b s - Again, assuming 
only free- free absorption is present, we have the condition 



would simply bounce off the clouds without being absorbed. 
However, the bounce itself is an impulse on the cloud that 
might have detectable effects. Thus, the radiation pressure 
of radio waves below the plasma frequency can squeeze IGM 
clouds, precisely because they reflect off them. 

In fact, the radiation pressure of radio waves can 
also squeeze t he IGM clouds if radio waves are effectively 
absorbed by (jCiotti fc Ostrikerl l2004h or scattered within 
them. Even without scattering, radio waves will be absorbed 
by free- free absorption for frequencies lower than v ca ,\(R) 
(eq. [37}. Scattering can only increase the coupling between 
radio waves and the cloud. If the naive scattering mean free 
paths in equation [2S1 ar e used, an IGM cloud is optically 
thick to scattering below a frequency 



(47) 



I conservatively assume that there is no scattering, so that 
the maximum frequency z^ nt at which a cloud is optically 
thick to either absorption or scattering simply equals u ca ,\(R) 
(eq.|37). 

There are two characteristic speeds for the collapse. The 
first is the sound speed of the IGM (c.f. eq. [7}, which is 



10/cT 



3mj 



1.7 x 10° cm s 



-lml/2 
1 i ■ 



(48) 



A relatively small pressure imbalance will cause the cloud to 
collapse, essentially quasi-statically and adiabatically, until 
its internal pressure is enough to resist the imbalance. If 
the sound crossing time is small compared to £igm ? then 
the pressure imbalance would have forced an IGM cloud 
to have a structure different than the one we observe. If 
R/c s < tiGM, the external sub-up radiation pressure P ra d 
cannot be much greater than the thermal pressure of the 
cloud P c m 2n e kT, or 



w ra d<2.0 x 10 18 ergs cm 3 S e Ta. 



(49) 



Ave,, < 



2nkT t abs 
tiGM tn 



< 1.6 x 10 10 L Mpc 



-3, ,2 r ,p5/2 



tiGM tft 



Gyr 10 Gyr 



(46) 



where again, T4 = T/(10 4 K). These bounds are naturally 
much weaker, since the IGM is transparent to radio waves 
at high enough frequency, and the radio waves keep most of 
their energy. 



5 RADIATION PRESSURE ON IGM CLOUDS 

No radio waves below the plasma frequency can traverse 
an IGM cloud. Since most of the IGM by volume is ex- 
tremely low density (5 < 0.1) (|Bi fc Davidsenl 119971 ). this 
means that most of the IGM could be filled with extremely 
low frequency waves (y < Hz) which cannot enter the IGM 
clouds that serve as our detectors in § [4] These radio waves 



8 Since we care about the absorption time and not the displace- 
ment from the source before absorption, the optical depth we care 
about is the true absorption optical depth, not the effective ab- 
sorption optical depth. 



When the external pressure is much greater than the 
(45) cloud pressure, however, the collapse can occur much more 
_.. quickly. Consider what would happen at the edge of a cloud 
immersed in a sub-i^ nt radiation bath with an extremely 
large pressure P ra d- The radiation pressure will rapidly ac- 
celerate the particles in the edge of the cloud until they 
exceed the local sound speed. A shock will form in the edge 
of the cloud. The shocking will continue until the local pres- 
sure balances the external pressure. Either the post-shocked 
gas will continue interacting with the sub-i/i nt radiation, in 
which case the exterior radiation pressure will continue to 
drive it inwards, or it will be transparent, in which case the 
radiation will now shock the next layer in until the entire 
cloud is shocked. 

Therefore, a shock will propagate inwards through the 
cloud, heating it until the internal pressure equals the ex- 
ternal pressure P ra d- The Rankine-Hugoniot conditions for 
a shock imply 



^ = 1 + Jl^ 
Pi 1 9 + 1 



{Ml - 1), 



(50) 



where Pi is the pre-shock pressure, P2 is the post-shock 
pressure - equal to P ra d in this case, Mi = v/c s ,i is the 
shock Mach number through the unshocked gas, and y g is 
5/3 for an ionized monatomic gas. We find that 
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Table 3. The IGM as a Free-Free Calorimeter 



Phase 


R 


v P 




^cal(etiGM) 6 


ve v Upper Bound c 




(Mpc) 


(Hz) 


(Hz) 


(kHz) 


(L Mpc- 3 ) 


DdCKgl U UI1U. IVjxlvl 


1 nno 

1UUU 


n 9 




u.uzu 




T T t*i t~*> v r i rincp T mroot 
V_JIlLLClU.t_.110C 1_jVCi 1CJ1 LB l 


I 


1 4 


7.5 


0.93 


5600 


Lya forest cloud 


1 


4.5 


72 


8.9 


56000 


Low density WHIM 


1 


10 


U 


1.7 


2.8 x 10 7 


High density WHIM 


1 


45 


780 


31 


5.6 x 10 8 


Cool WHIM 


1 


9.0 


870 


42 


1.5 x 10 5 


Low density cluster 11 


1 


90 


110 


1.5 


2.3 x 10 12 


High density cluster 1 * 


1 


280 


1100 


14 


2.3 x 10 13 


Strong Mg II absorber 


0.001 


900 


2.4 x 10 5 


17000 


6.7 x 10 11 


Dense weak Mg II absorber 


10~ 5 


900 


26000 


3100 


2.3 x 10 13 


Weak Mg II/C IV absorber 


0.001 


280 


26000 


3100 


2.3 x 10 10 


C IV/O VI Galaxy Halo 


0.1 


28 


2700 


140 


1.4 x 10 7 



a Frequency at which the mean free absorption path is the width of the IGM phase, the minimum possible 
value for i/ ca ] (if only free-free absorption is present). 

h Frequency at which the free- free absorption time-scale is tiGMi the maximum possible value for v ca \. 
c Upper bound on the luminosity density within each cloud for v < v c!l \. Assumes that the luminosity is 
distributed continuously throughout the cloud, so that no part of the cloud is shielded, and that there are 
no other sources of absorption. 

d Synchrotron absorption may exceed free-free absorption, making the luminosity density bound invalid for 
v > vr- 



Pad = \ P iv 2 - ^Pt (51) 

In the case when P ra d = Pi, we simply get v = c s ,i as 
expected. Since the IGM clouds we do see have not been 
destroyed yet, we can suppose that v < R/tiGM- Using the 
fact that u ra d = 3P ra d, we find that 

"rad < [9.0 x 10~ 17 <5 e PMp C tro 2 - 2.6 x 10~ 19 <5 e T 4 ] ergs cm~ 3 (52) 

where tio = £igm/10 Gyr, Pm P c = P/Mpc, and T4 = 
T/(10 4 K) for the initial cloud. 

If the post-shock region is radiating efficiently, the shock 
may instead be isothermal, with the final temperature equal 
to the initial temperature. In this case 

^ = Mi, (53) 
r\ 

where P3 is the pressure of the radiatively cooled region. If 
we now assume that P3 = P ra d, we find that 



^rad 



< 7.2 x 10 17 <S e PM pc f 10 2 ergs cm 3 . 



(54) 



In all considered cases, the limits on ti ra d from the isothermal 
case are not as conservative as those from the non-radiating 
case. 

Equation[S5]holds only if Mi > 1, because of the Second 
Law of Thermodynamics; in the case when P/tiGM < c s ,i, 
equation [49] is appropriate since pressure imbalances will be 
propagated at least at the sound speed. However, since £igm 
was chosen to be at most t SOU nd in § 12.21 eq.[52]is always the 
appropriate choice. 

Note that the radiation pressure bound is bolometric for 
all electromagnetic radiation below v- m %. All incident radio 
waves with frequencies below this frequency will be reflected 
off, scattered by, or absorbed by the cloud, and will there- 
fore squeeze the cloud. Only if the wavelength of the incident 
waves are greater than the size of the cloud will the radi- 
ation fail to interact, and the background plasma density 
with a low enough plasma frequency for a radio wave with 
wavelengths of kpc to Mpc to exist is inconceivably low. 



6 SYNCHROTRON HEATING IN CLUSTERS 

Cluster gas is very hot and dense, so the limits I derived 
from free-free absorption (§ 14. 1[) and radiation pressure (§[!} 
for the sub-MHz radio background in galaxy clusters are 
weak. However, there is another gas phase in galaxy clusters 
that can serve as a sub-MHz detector: the relativistic phase. 
Some clusters are filled with high energy cosmic ray electrons 
and possibly positrons (CRs), which emit sync hrotron radio 
emission observed at MHz to GHz frequencies |Ferrari et al.l 
120081 ). These electrons also Inverse Compto n scatter CMB 
photo ns to produce hard X-ray emission (Rcphacli ct al. 
2008). Synchrotron cooling and IC cooling off the CMB 
should dominate the cooling processes at high energies. In- 
verse Compton losses off starlight produced by the cluster 
are relatively weak; given a typical cluster luminosity of L ~ 
10 46 ergs s _1 and a radius of R ~ 1 Mpc, the typical starlight 
energy density is = L/(A-kR 2 c) Ri 3 x 10~ 15 erg cm -3 . 
Indeed, cluster radio emission is observed to fall off steeply, 
as would be expected by synchrotron and IC cooling. 

However, a large radio background can actually ac- 
celerate particles, especially at lower energies, through 
synchrotron absorption. This effect can overcome Inverse 
Compton and synchrotron cooling in ce rtain circumstances 
l|Ghisellini. Guilbert. fc Svenssonl [1988). The synchrotron 
heating rate is 



bhe 







-2 E dE 



N(E) 



E 2 



N{E)~ 



k;P(v,E)dv(tt) 



Synchrotron heating will create a rising CR elec- 
tron/positron spectrum at low energies that peaks at some 
frequency depending on the IC and synchrotron cooling 
time-scales. Therefore, the CR electrons and positrons in 
a galaxy cluster also serve as a radio detector. The very fact 
that no such peak is seen in the radio spectra of clusters 
places limits on the sub-MHz radio background in these re- 
gions. 

There is a competition between heating and cooling 
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Table 4. Radiation Pressure Constraints on sub- 



backgrounds. 



Phase 


fcal(-R) 




pb 

therm 




Cs 


Mrad(^ < Vint) Upper Bound c 




(Hz) 


(kHz) 


(erg cm -3 ) 


(km s _1 ) 


(km s" 1 ) 


(erg cm -3 ) 


~r~> i 1 T^-| Ti jr 

Background 1GM 


15 


11 


"i a in — 21 

1.4 X 10 


f\ o c\c\c\ 

98000 


17 


1.8 X 10 


Underdense Lya forest 


7.5 


52 


f c\ v , "i r\ — -20 

6.9 x 10 


160 


17 


o -i in — 1 Q 

2.1 X 10 


Lya forest cloud 


72 


180 


6.9 x 10~ 19 


160 


17 


2.1 x lO" 18 


Low density WHIM 


11 


540 


3.5 x 10~ 16 


160 


170 


1.0 x 10~ 15 


High density WHIM 


780 


2800 


6.9 x 10" 15 


160 


170 


2.1 x 10~ 14 


Cool WHIM 


870 


480 


2.8 x 10" 18 


11 


17 


4.4 x lO" 16 


Low density cluster 


110 


5900 


2.8 x 10~ 12 


1600 


1700 


8.3 x 10~ 12 


High density cluster 


1100 


2.1 x 10 4 


2.8 x lO" 11 


1600 


1700 


8.3 x lO" 11 


Strong Mg II absorber 


2.4 x 10 5 


3.9 x 10 4 


2.8 x 10~ 14 


49 


17 


8.8 x 10~ 13 


Dense weak Mg II absorber 


26000 


2.5 x 10 4 


2.8 x 10" 14 


16 


17 


8.3 x 10~ 14 


Weak Mg II/C IV absorber 


26000 


1.1 x 10 4 


2.8 x 10" 15 


16 


17 


8.3 x 10~ 15 


C IV/O VI Galaxy Halo 


2700 


1400 


2.8 x 10~ 17 


98 


17 


3.6 x 10~ 15 



a Frequency at which the cloud becomes optically thick to scattering, if eg. 1201 describes the scattering mean free path and A n = 1 
and £q = R. 

b Pressure of the cloud; if R/tiQM < c s , this is the limit on the external sub-i/p radiation pressure. 
c Note that « ra[ j = 3P ra( j. 



(including Inverse Compton cooling off the low frequency 
background itself) at every CR energy. In order for syn- 
chrotron heating to be effective, the absorbed radiation fre- 
quency i/o must be greater than the Razin frequency vr 
(eq. [3]), and the synchrotron heating time-scale must be 
less than the cooling time-scale from all other processes: 

ihcat < [£ synch + ^IC,CMB + ^IC, sub-MHz + *IGm1 ' ^ nC S y n ~ 

chrotron cooling time-scale of electrons/positrons observed 
in synchrotron radio at vc = /^c,GHzGHz is 



t S ynch = 1.59 x 10 yr v y 



-1/2 



C.GHz^/jG 



3/2 



(56) 



The Inverse Compton cooling time-scale from low frequency 
photons at frequency near uo is 



tic = 6.32 x 1(T 5 yr ^ c ^Xg [ A ^M]~\ 



(57) 



where Avou u (vo) is in cgs (ergs cm ). Additional Inverse 
Compton cooling comes from the CMB: 



(58) 



tic.CMB = 1.50 x 10 8 yr v c %hz B Ig- 

The synchrotron heating time-scale off photons near fre- 
quency vr < Vo -C v, if the electron spectrum goes afl E~ 3 , 



t hca t = 1.93 x 10~ 7 yr u c %^' 3 Kz B^v o u v {u )\-\ (59) 

As seen in Figure [21 if the sub-MHz radio background is 
the maximum allowed by the free-free absorption bounds, 
synchrotron heating is by far the most important process at 
low energies for a galaxy cluster. The fact that we do not 
see a turnover in the radio synchrotron spectrum of galaxy 
clusters at ~ 10 GHz implies that the sub-MHz radio back- 
ground in clusters must be far smaller than that allowed 
by the free-free absorption bounds. Using these time-scales 
gives us the following condition for synchrotron heating to 
be too weak to alter the spectrum: 



9 There is some residual dependence on the electron spectrum 
shape, because the effectiveness of synchrotron absorption is af- 
fected by stimulated synchrotron emission. Using an E~ 2 or E~ 4 
CR electron/positron spectrum will not alter these estimates 
much. 



Az;o^(fo) < 10 16 ergscm 



l.2Bt G + 13 + 1.9S; G £ 



1/2,-1, 1/2 
nG L 9 C,GHz 



-11/6 -5/3 fi 5/2 
^C.GHz^O.kHz^G 



0.000305 



(60) 



where t 9 = iiGM/Gyr. 

In practice, the synchrotron heating bound appears 
three-sided, as seen in Figure [3] Synchrotron heating by low 
frequency radio waves with uq < vr is not effective because 
of the Razin cutoff. At high CR electron/positron energies, 
Inverse Compton cooling off the low frequency radio back- 
ground becomes more effective than synchrotron heating 
from those same photons. Then synchrotron heating ceases 
to be important. The target photon frequency Vq where this 
happens for electrons/positrons observed at ^c.GHzGHz is: 



oil tr -11/10 D 3/2 

vo = 32kHz v c ( J^ BJ G . 



(61) 



This is the high-frequency cutoff for the excluded region 
visible in Figure [3] Finally, for intermediate frequencies vq, 
the synchrotron heating must primarily compete with IC 
cooling off the CMB, and to a lesser extent, synchrotron 
cooling and adiabatic losses. We see from equation 1601 that 
the bounds on AvoUv(iso) in this regime go as u^/ 3 , as the 
shown for bottom of the excluded region in Figure [3] 

For a galaxy cluster with microGauss magnetic fields, 
the upper limit on radio emission above the Razin cutoff is 
very small from the lack of a spectral peak at GHz: 



<10~ 



ergs cm 



11/6 5/3 



x C,GHz' 



0,kHz-^fjG 



5/2 



(62) 



The synchrotron spectrum of the Coma cluster has al- 
ready been observed down to 30 MHz, and the limits on 
its sub-MHz energy density are strong. LOFAR, which can 
observe all the way down to 15 MHz, can place incred- 
ibly strong bounds on kHz to MHz radio emission. For 
example, synchrotron heating from 100 kHz radio would 
alter the radio spectrum at these frequencies as long as 
vqu v (vq) > 1 X 10 -15 ergs cm -3 B~q 2 . This is roughly com- 
parable to the starlight energy density of a galaxy cluster; 
recall that, unlike starlight, sub-MHz radio emission might 
be trapped for a long time in clusters by plasma scattering, 
so that a very low luminosity can accumulate for a long time 
to give a large energy density. 

The synchrotron absorption argument should also ap- 
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Figure 2. Cooling and heating times for CR in a galaxy clus- 
ter with a magnetic field strength of 0.5 /iG and a gas density 
of 10 -3 cm -3 . For the synchrotron heating times, I assume each 
AvqUv(vo) is the free-free absorption upper limit on a homoge- 
neous background filling the cluster gas (§ 14,111 . for uo = 100 kHz 
and 1 MHz. Although synchrotron heating by 100 kHz photons is 
more effective than by MHz photons at fixed Ai/qu^^o), the up- 
per limit on the 100 kHz background is much more severe, so the 
synchrotron heating allowed by the free-free absorption bound is 
weaker. 



ply wherever there are strong enough magnetic fields and 
observable cosmic ray electrons. For example, it should also 
apply to structure formation shocks, which are believed to 
accelerate cosmic rays and produce radio emission. 

In addition to syn chrotron absorpti on, there should be 
transition absorption ([Fleishm an 1989). In analogy with 
synchrotron heating, we would expect a transition heating 
effect, where low energy particles are heated by transition 
absorption. I calculated the transition heating time-scale for 
an E~ 3 CR electron/positron spectrum, using equation 1551 
Unlike the synchrotron reheating times, transition absorp- 
tion is inefficient in clusters: 

theat = 4.5xl0 6 yr B~QU C ,GH^l^ Hz n^£l / ^ [pc A~ 1 [u u v (i'„)r L 

Therefore transition absorption provides no strong con- 
straints on the sub-MHz background in clusters. 



7 INVERSE COMPTON BOUNDS IN 
CLUSTERS 

The CR population of electrons in galaxy clusters should 
also Inverse Compton scatter any sub-MHz radio emission 
in the cluster to much higher frequencies, such as MHz to 
GHz radio. Therefore the MHz to GHz radio emission of the 
cluster limits the amount of sub-MHz radio of the cluster. 

The Inverse Compton emissivity of photons of energy 
E 1 from an electron/positron spectrum N e (E e ) scattering a 
photon field with number density n D h(e) at original photon 
energy e is given by (e.g., Schlickeiscr 2002) as 



Figure 3. Limits on the homogeneous sub-MHz energy densities 
in a cluster with n = 10 -3 cm -3 , B = 0.5 fJ,G, T = 10 8 K, 
and tiGM = 0-6 Gyr. The dashed line [light grey shading) is the 
(optically thin) free- free absorption limit (cq. [31}. Solid lines are 
synchrotron heating limits, assuming that there is no spectral 
downturn at synchrotron frequencies of 10 MHz to GHz. Syn- 
chrotron reheating limits are generally much stronger than the 
free-free absorption limits, but only apply above the Razin cutoff 
{vr) and require a low frequency radio detection of the cluster. 
The Inverse Compton limits (dotted) for a cluster with a Coma 
cluster-like CR electron spectrum are the strongest constraints 
below the Razin cutoff. 



da(E 1 , e, E e 



■cN e (E e )dE e , (64) 



<2 7 (£ 7 )ic = / n ph (e)de 

in photons per unit volume per unit time per unit energy. 
The cross section in this formula is 



da S(7t , 9, 
^ = 4^ (29ln « +1 + *- 29) 



(65) 



and q « i? 7 /(47 2 e) in the Thomson limit. The Inverse 
Compton emissivity at frequency v\ — E~,/h is converted 
into a specific luminosity as L v — h 2 v\VQ 1 (E 1 ). The con- 
tion that L v < L„, bs at frequency v\ then becomes 



'(yo) 



4L 

3&tcV v\ 



(2glnq + 1 + q 



2,dN _ 



In Figure [3l I show the Inverse Compton bounds on the 
sub-MHz emission in the Coma Cluster (dotted), using the 
CR electron spectrum derived from the synchrotron MHz to 
GHz radio spectrum (§ 12. 3[) . Near the plasma frequency, the 
IC bound is 8 x 10 -15 erg cm -3 , which is somewhat greater 
than the energy density in starlight in a cluster. At higher 
frequencies, the IC bound grows as j~ 2 . An electron that 
emits 30 MHz synchrotron radiation (7 ~ 1000) can boost 
kHz photons to 4 GHz photons, the highest frequency radio 
emission seen from the Coma cluster. When va > kHz, the 
integral over dN/d'y is taken to be constant, since we do not 
know the spectrum at lower energies, and every CR elec- 
tron/positron observed in synchrotron is energetic enough 
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to boost the v>o frequency photons to observable photons 
of frequency V\ E3 Thus the bound on u v (vq) goes as Vq. If 
the low energy electron spectrum is ever known better, that 
would increase the strength of the bounds at high v§. 

We see that the IC bounds are stronger than the free- 
free absorption bound (dashed), which were weak because 
cluster plasma is dense and hot. The Inverse Compton con- 
straints are not as strong as the synchrotron heating con- 
straints (solid) in a cluster like Coma, where data below 
100 MHz is available. Unlike the synchrotron heating con- 
straints, the Inverse Compton bound extends below the 
Razin cutoff. 

LOFAR will be able to measure the radio spectrum of 
several clusters down to ~ 15 MHz. As with the synchrotron 
heating bounds, these radio observations will allow us to 
set strong limits on the sub-MHz radio emission in several 
clusters by the Inverse Compton argument. 



8 POTENTIAL LIMITS FROM ULTRA HIGH 
ENERGY 7-RAYS 

If two photons each have an energy greater than m e c 2 in 
their centre of mass frame, the pair of photons can convert 
into an electron-positron pair. For photons of energy Ei and 
E2 in the observer frame, the pair-production condition is 
that E1E2 > m 2 c 4 . Pair production is expected to be an im- 
portant source of opacity for PeV pho tons as they interact 
with the CMB (Mo skalenko et alfe OOG). A large radio back- 
ground can also serve as a source of opacity for ultra high 
energy (UHE) photons. For a MHz radio photon, however, 
the target photons must have energy 2 x 10 20 eV, and for 
kHz radio, the target photons must be at least 2 x 10 23 eV. 

No UHE photons have ever been detected, but if extra- 
galactic photons of high enough energy were ever detected, 
they would provide potent constraints on the sub-MHz ra- 
dio background. The absorption length scale for pair produc- 
tion is A 77 ~ (n-yOT)^ 1 ■ The path length is 1 Mpc only when 
n 7 ~ 0.5 cm -3 or it « 3.2x 10 _24 ^hz erg cm -3 . For compar- 
ison, the 5 = 1 Lya forest free-free absorption bound implies 
the sub-MHz radio photon density is n 7 ~ 3300 fkHzCm -3 
even in the optically thin, homogeneous case. Thus detec- 
tions of extragalactic photons with energies above 10 20 eV 
would increase the bounds on the extragalactic sub-MHz 
radio background by several orders of magnitude. 

Pair production limits would have the advantage of 
being relatively more powerful at higher frequency, as op- 
posed to the free-free absorption bound which is more 
powerful at lower frequency. It is more likely that UHE 
photons with lower energy will be detected if their spec- 
trum is falling, which in turn requires higher energy ra- 
dio photons for pair production to happen. Indeed, there 
probably is some flux of photons of energy 10 20-21 eV 
produced by the Gre i den-Zatsepin-Kuzmi n process (e.g., 
IWdowczvk et al.lll972l ; iGelmini et ai1l200Sl ). while there is 
no obvious mechanism for making photons of much higher 
energy than that. Also, the pair production path length 
scales inversely with the number density of photons, or 

10 The (2olno+l+q — 2a 2 ) in the integral of eq, l66l is of order ~ 1 
for vq >kHz for all CR electrons/positrons observed at 30 MHz < 
vi < 4~GHz. 



A 77 oc (/S.vu„)~ 1 v, whereas the free-free absorption path 
length scales as Xg oc (£±vu v )~ x v 2 ' . 

Unfortunately, there is a limit to these arguments. 
The Galaxy becomes completely opaque to photons with 
energy of 10 24 eV and higher, because single photons can 
pair produce of f the Galact ic magnetic field (|Steckerll200l ; 
see also § 3 of lErberlll966h . Therefore, we can expect no 
pair-production bounds on the extragalactic radio back- 
ground below about 250 Hz. Any such bounds also assume 
Lorentz invariance holds a t these high energies; some theo- 
ries do not predict this (see lAmelino-Camelia fe Piranll200ll ; 
iGalaverni fc Siglll2008l ). 



9 CONCLUSION 

The extragalactic sub-MHz background is invisible to direct 
observation from Earth, but we can still detect its effects on 
intergalactic matter. I have placed new limits on the magni- 
tude of the sub-MHz radio background, using various IGM 
phases and clusters as radio detectors. Figure [5] (full legend 
in Table [5]) summarises the bounds on the radio background 
from the IGM thermal state (§ 14.11 red) and the radiation 
pressure exerted on the IGM (§0 brown). A sub-MHz back- 
ground with an energy density as large as the CMB at any 
frequency is easily ruled out, and energy densities compara- 
ble to the cosmic starlight backgrounds are also not allowed 
at almost all frequencies. 

Low frequency radio waves can heat the IGM through 
free-free absorption. Observations of the IGM thermal state 
constrain the amount of heating from the extragalactic sub- 
MHz radio background. Free-free absorption bounds (§ [4TTJ 
are potentially the strongest of all of the limits, but are 
highly model dependent at low frequencies. If we are con- 
sidering a radio bath that pervades the entire cloud evenly 
(dashed red lines), then the entire cloud is heated up and the 
bounds are very strong at low frequency in terms of energy 
density. In fact, the energy density constraints within the 
Lyman-a forest just above its plasma frequency would be the 
strongest of any photon energy, as seen in Figure [4] If a back- 
ground was simply incident on the IGM clouds, then when 
the cloud becomes optically thick, the outside of the cloud 
will be heated but the interior will not. In this case, if there is 
no scattering of radio waves within the cloud (solid red lines, 
light pink shading), the clouds usually remain optically thin 
down to sub-kHz frequency, and the energy density bounds 
remain strong. However, scattering will increase the effec- 
tive absorption optical depth. Naively applying the scatter- 
ing mean free path in equation [20] considerably weakens the 
energy density bounds (dotted red lines, darker pink shad- 
ing). However, the approximations in equa tion | 20 | mav break 
down at high scattering optical depth (|2.3llCohen fc Cronvnl 
1974). I have also not considered the opacity of any remain- 
ing shell of 'evaporated' material around an optically thick 
cloud. In order to set firm limits from free-free absorption, 
we need to understand the radiative transfer of sub-MHz 
radio waves through the IGM better. 

I have also set an upper limit on the luminosity density 
within each IGM phase from free-free absorption (§ 14. 2[) . 
assuming the radio waves fill the IGM cloud evenly. From 
the existence and temperature of the Lyman-a forest, I in- 
fer that at frequencies of 5 — 100 Hz, these clouds have a 
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log 10 [E (eV)] 
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log 10 [v (Hz)] 

Figure 4. The z = cosmic backgrounds for the electromagnetic spectrum. The cosmic plasma frequency (assuming 8 e ~ 0.002) makes 
propagation impossible a t the lowest frequencies (dark grey); we expect there to be no radio background below this cutoff. The flu bound 
(Zentner & Walker 2002) is in grey at top (§01. Free-free absorption bounds (§|4j from Lya clouds of 8 = 0.1 and 1.0 and from a weak Mg 
II/C IV absorber are shown in pink/red, assuming J„ is a bump at v. The solid red bounds are for the case when the radiation fills the 
cloud evenly. The shaded regions are bounds on an incident external radiation field; lighter shading (dashed boundaries) for no scattering 
in the cloud, while darker shading (dotted boundaries) when eq. l20l describe the scattering. Radiation pressure can crush a cloud (§[5j; the 
bounds on cloud crushing are shown in brown. The solid line and shading assumes there is no scattering, while the dotted line assumes 
that eo. 1201 describe the scattering. The radiation pressure bounds are bolometric below Vi nt . I also plot naive upper expectations on 
the radio background expected from several sources (S 12.111 . not accounting for IGM absorption: the maximum synchrotron brightness 
temperature (Tj £S 10 12 K), pulsars, and gravitational wave conversion. See Table [5] for a full legend with references. 
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Table 5. Legend for Figure [1] 


Wavelength band 


Method/Instrument 


Reference 


Symbol 


All 


O 


ZiCntncr & Walker ( zUUz ) 


light grey shading 


bub-JVlriz radio 


IGM pressure 


1 nis work lop 


brown shading 




IGM thermal state 


inis work ( q 14. ip 


pink shading 


Radio 


Theoretical prediction 


Frotncroc & Bicrmann (lyyo) 


dotted black lines 




A T~i /"< A TNTTIO 

ARCADEi 


lixscn ct al. (2009) 


open triangles 


IR to UV 


Theoretical prediction 


Iranccslimi, Rodigincro, & Vaccari (2008) 


solid black line 


Infrared 


T?T"D A Q 
£ IKAb 


r lxscn ct al. ( lyyo ) 


orange dotted line 




TiT A c~n 

BLAb 1 


Marsdcn ct al. (2009) 


orange circles 




DIRBE 


Tier ■ i j_ / i~\r\r\ a \ 

Wright (2004) 


orange pentagons 




Spitzer 


Dole ct al. (2006) 


green lower limits 




Spitzer 


Papovich et al. (2004) 


orange triangle 




Spitzer 


Savage & Oliver (2005) 


green squares 


Optical 


Galaxy counts with Hubble 


IVladau & Fozzctti (zUUUJ 


filled triangles 


U V 


PAT T?"Y" 


yS.ll ct al. ( ZUUO ) 


violet triangles 


ri/U V 


LyQ forest ionization 


briuii ct ai. (lyyyj 


open circle 


X-rays 


XMM-Newton (Lockman hole) 


Worslcv ct al. (zUUoJ 


/ _ i 

grey o-stars 




Chandra 


TT" 1 TIT 1 "j_ 1 /ririM^N 

fiickox & Markcvitch (2006) 


5-stars 




Swift 


Moretti ct al. (2009) 


dotted blue line 




Swift 


Aicllo ct al. (2008) 


solid blue line 






Kcvmvtscv ct al. (zUUoJ 


solid green line 




HEAOl 


Kmzcr et al. (1997) 


grey solid line 


MeV 7- rays 


SMM 


Watanabe et al. (2000) 


solid cyan line 




COMPTEL 


Wcidcnspointner et al. (2000) 


crosses 


GeV 7- rays 


EGRET 


Strong. Moskalenko. & Reimer (2004) 


open squares 


> TeV 7-rays 


GeV background 


Coppi & Aharonian (1997) 


violet shading 


TeV 7-rays 


HESS 


Aharonian et al. (2009) 


solid black line 




HESS 


Aharonian et al. (2008) 


solid violet line 




GRAPES-3 


Havashi & Grapes Collaboration (2003) 


blue arrows 


PeV 7-rays 


CASA-MIA 


Chantell et al. (1997) 


grey arrows 


EeV 7-rays 


Auger 


Abraham et al. (20091 


black arrows 




Auger 


Abraham et al. (2008a) 


black arrows 



The PeV 7-ray limits do not correct for pair-production absorption from the CMB. 



maximum emissivity of 6 x 10 4 Lq Mpc -3 . Again, a better 
understanding of the radiative transfer and the scattering 
in particular of sub-MHz radio waves is needed to set more 
firm limits on the luminosity density of the Universe at low 
frequencies. 

At the lowest frequencies, there is a window in the free- 
free absorption constraints. This is because voids are ex- 
tremely underdense, with a very low plasma frequency. Sub- 
Hz radio could be generated in the voids and would simply 
reflect off the more condensed structures that would other- 
wise be heated by them. 

Free-free absorption constraints are weak because scat- 
tering or absorption of the radio waves shields the interior of 
each IGM cloud. At the very lowest frequencies, below the 
plasma frequency of the cloud, the radio waves simply re- 
flect off it. There are none the less constraints even at these 
lowest frequencies, because the reflection, scattering, or ab- 
sorption of these waves squeezes IGM clouds (§[5]). These 
radiation pressure bounds (Figure U brown) are somewhat 
weak but are still strong enough to rule out an average sub- 
Hz background as large as the starlight backgrounds. Un- 
like the free-free absorption bounds, the maximum u v are 
not model dependent. However, the frequency range over 
which they are applicable also depends on the radio scatter- 
ing properties of the IGM. At the very least, the radiation 



pressure bounds apply until the cloud is optically thin to 
free-free absorption. 

Galaxy clusters have hot and dense gas, which makes 
them poor free-free absorption detectors for any sub-MHz 
radio waves within them. However, they also contain mag- 
netic fields and cosmic rays, which can interact with sub- 
MHz radio waves in additional ways and provide additional 
limits (Figure [3]). Low frequency radio waves above the 
Razin frequency (~ 20 kHz) can actually heat CR elec- 
trons/positrons, and would create a peak into the observed 
MHz to GHz synchrotron radio spectra of clusters (§[6]). The 
lack of such a peak rules out sub-MHz radio backgrounds as 
small as ~ 10 _15 ergs cm -3 in the Coma cluster, compara- 
ble to the energy density in starlight. CR electrons/positrons 
can also Inverse Compton scatter low frequency radio waves 
to observable MHz to GHz frequencies (§0. The observed 
radio spectrum again constrains the sub-MHz radio back- 
ground in Coma to be as small as ~ 10~ 15 ergs cm -3 for 
!^<kHz. 

Some relatively weak statements can be made about 
whether the backgrounds described in § I2.1l exist. The most 
exotic sources of sub-MHz radio waves are constrained. The 
radiation pressure bound from weak Mg II/C IV absorbers 
are sufficiently strong to exclude a radio background as large 
as the LIGO upper limits on a 100 Hz stochastic gravi- 
tational wave background. A background of pulsar waves, 
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if they somehow escaped into the IGM and did not suffer 
absorption, would be weaker still; the free-free absorption 
bounds without scattering are strong enough to rule out 
such backgrounds from all young pulsars and all MSPs in 
Galactic discs. These maximum estimates of the sub-MHz 
radio background are probably unrealistic anyway (see the 
discussion in § I2.1.3[) . Unfortunately, a more realistic syn- 
chrotron background also seems to be out of reach by the 
free-free absorption bounds; a Tj < 10 12 K background is 
ruled out only for i/>MHz, which can already be directly 
observed. However, if the scattering properties of the IGM 
are similar to those in the ISM, such that eq. I20l holds. then 
the radiation pressure bounds from the Lyman-a forest will 
rule out maximal synchrotron backgrounds at ~ 100 kHz. 

There are several ways to reduce the uncertainties in 
these bounds. Knowledge of the low frequency scattering 
properties of the IGM is essential for the free-free absorp- 
tion bounds, which would otherwise be strong (solid red lines 
in Figure |4}. This knowledge can also help us determine 
the frequency range the radiation pressure bounds apply 
over. Strong scattering will weaken the free-free absorption 
bounds, because radiation cannot diffuse deep into the cloud 
and heat it; but strong scattering strengthens the radiation 
pressure bounds, because radiation can then efficiently cou- 
ple with the cloud exterior and squeeze it. Low frequency ra- 
dio observations of galaxy clusters and other environments 
with CRs are especially useful for the Inverse Compton and 
synchrotron heating bounds. Observations of extragalactic 
ultra high energy photons, if they exist, would set extremely 
strong constraints on the extragalactic sub-MHz background 
(§[8]), especially at high frequencies where the other bounds 
are weakest. 

It would be a simple matter to apply similar constraints 
to sub-MHz emission within the Galaxy itself, which will be 
done in a future paper. Not only are the thermal proper- 
ties of the interstellar medium relatively well characterised, 
but the Galaxy has a well known CR electron spectrum and 
magnetic field. Therefore, we could apply synchrotron heat- 
ing and IC upscattering arguments, which are not dependent 
on the scattering of low frequency radio waves. This would 
allow us to probe distant regions of the Galaxy that are 
not visible at low radio frequencies because of free-free ab- 
sorption, such as the Galactic Centre, which is obscured by 
free-free absorptio n at frequencies as high as 330 MHz (e.g., 
IPedlar et al1ll989h . 

Although the bounds in this paper may not strongly 
constrain expected sources like a synchrotron background, 
the extragalactic sub-MHz sky is not completely unknow- 
able. Instead of disregarding the low frequency emission of 
radio sources, it is possible to consider the effects of the 
emission on their surroundings. These effects may prove to 
be important to our understanding of the regions around 
sub-MHz sources. Even if the extragalactic sub-MHz sky is 
forever invisible to us directly, its presence can still be seen. 
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